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^ ■ Abstract 

. We perform statistical analysis of the single-vehicle data measured on the Dutch freeway A9 and dis- 

cussed in 0. Using tools originating from the Random Matrix Theory we show that the significant changes 
Q , in the statistics of the traffic data can be explained applying equilibrium statistical physics of interacting 

' particles. 
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^ ■ The detailed understanding of the processes acting in the traffic systems is one of the most essential parts of 
7— ( . the traffic research. The basic knowledge of the vehicular interactions can be found by means of the statistical 
analysis of the single-vehicle data. As reported in Ref. Q, 0, and 0] the microscopical traffic structure can 
' be described with the help of a repulsive potential describing the mutual interaction between successive cars 
""^ ■ in the chain of vehicles. Especially, the probability density Pp(r) for the distance r of the two subsequent cars 
O | (clearance distribution) can be described with the help of an one-dimensional gas having an inverse temperature 
c/2 ' (3 and interacting by a repulsive potential V(r) = r" 1 (as discussed in Ref. ||5l and 121). Such a potential leads 
to a clearance distribution 

Pp(r)=Ae-%- Br , (1) 
J> . where the constants A = A(f3), B = B(j3) fix up the proper normalization Pp(r) dr = 1 and scaling 
^ ! Jo°° r Pp( r ) dr = 1. This distribution is in an excellent agreement with the clearance distribution of real -road 
$— i ' data whereas the inverse temperature (3 is related to the traffic density q. 
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Another way to seek for the interaction between cars within the highway data is to investigate the traffic 
flow fluctuations. One possibility is to use the so-called time-gap variance At considered in paper [23 and 
defined as follows. Let {U : i = 1 . . . Q} be the data set of time intervals between subsequent cars passing a 
fixed point on the highway. Using it one can calculate the moving average 

k+N-l 

T * N) = N £ * (k = l...Q-N+l) 
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Figure 1: The time-gap variance At(N) as a function of the sampling size N (in 
log-log scale). Plus signs and stars represent the variance of average time-gaps for free 
and congested flows, respectively. 



of the time intervals produced by the N + 1 successive vehicles (i.e. N gaps) as well as the total average 
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The time-gap variance A^ is denned by the variance of the sample-averaged time intervals as a function of the 
sampling size N, 
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where k runs over all possible samples of N + 1 successive cars. For time intervals ij being statistically inde- 
pendent the law of large numbers gives At{N) ocl/N. 

A statistical analysis of the data set recorded on the Dutch freeway A9 and published in Ref. |2]| leads, 
however, to different results - see the Figure 1. For the free traffic flow (g < 15 veh/km/lane) one observes 
indeed the expected behavior At(N) oc 1/N. More interesting behavior, nevertheless, is detected for higher 
densities (g > 35 veh/km/lane). Here Nishinari, Treiber, and Helbing (in Ref. J2)) have empirically found a 
power law dependence 

A r (AT) oc iV 7 

with an exponent 7 « —2/3, which can be explained as a manifestation of correlations between the queued 
vehicles in a congested traffic flow. 



There is, however, one substantial drawback of this description. The time-gap variance was introduced ad 
hoc and hardly anything is known about its exact mathematical properties in the case of interacting vehicles. It 
is therefore appropriate to look for an alternative that is mathematically well understood. A natural candidate 
is the number variance A n (L) that was originally introduced for describing the statistics of eigenvalues in 
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the Random Matrix Theory. It reproduces also the variances in the particle positions of certain class of one- 
dimensional interacting gases (for example Dyson gas in Ref. if!!). 

Consider a set {rj : i = 1 . . . Q} of distances (i.e. clearances in the traffic terminology) between each pair 
of cars moving in the same lane. We suppose that the mean distance taken over the complete set is re-scaled to 
one, i.e. 

Q 

i=l 

Dividing the interval [0, Q] into subintervals [{k — 1)L, kL] of a length L and denoting by nk{L) the number 
of cars in the fcth subinterval, the average value n(L) taken over all possible subintervals is 

1 IQ/ L \ 

= Wl\ £ nk{L) = L ' 

where the integer part [Q/L\ stands for the number of all subintervals [(k — 1)L, kL] included in the interval 
[0, Q]. Number variance A n (L) is then defined as 

, IQ/H 

- E ("»<*> - 1) 2 

and represents the statistical variance of the number of vehicles moving at the same time inside a fixed part of 
the road of a length L. The mathematical properties of the number variance are well understood. For indepen- 
dent events one gets A n (L) = L. Applying it to the highway data in the low density regime (free traffic) one 
obtains however A n (L) as 5L/6 (not plotted). The small deviation from the behavior A n (L) = L is induced 
by the weak (but still nonzero) interaction among the cars. 

The situation becomes more tricky when a congested traffic is investigated. The touchy point is that behav- 
ior of the number variance is sensitive to the temperature of the underlying gas - or in the terminology of the 
Random Matrix Theory - to the universality class of the random matrix ensemble. To use the known mathemat- 
ical results one has not to mix together states with different densities - a procedure known as data unfolding 
in the Random Matrix Theory. For the transportation this means than one cannot mix together traffic states 
with different traffic densities and hence with different vigilance of the drivers. So we will perform a separate 
analysis of the data-samples lying within short density intervals to prevent so the undesirable mixing of the 
different states. 

We divide the region of the measured densities q G [0, 85 veh/km/lane] into eighty five equidistant subin- 
tervals and analyze the data from each one of them separately. The number variance A n (L) evaluated with the 
data in a fixed density interval has a characteristic linear tail (see Fig. 2) that is well known from the Random 
Matrix Theory. Similarly, such a behavior was found in models of one-dimensional thermodynamical gases 
with the nearest-neighbor repulsion among the particles (see Ref. (6)). We remind that for the case where the 
interaction is not restricted to the nearest neighbors but includes all particles the number variance has typically 
a logarithmical tail - see 0D|. So the linear tail of A n (L) supports the view that in the traffic stream the interac- 
tions are restricted to the few nearest cars only. The slope of the linear tail of A n (L) decreases with the traffic 
density (see the top subplot in the Fig. 3). It is a consequence of the increasing alertness of the drivers and hence 
of the increasing coupling between the neighboring cars in the dense traffic flows. 
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Figure 2: Number variance A n (L) as a function of length L. Plus signs, diamonds, 
stars, circles, crosses, and squares represent the number variance of real-road data in 
selected density regions (see legend for details). The curves show the linear approxi- 
mations of the individual data. Their slopes were carefully analyzed and consecutively 
visualized on the Fig. 3 (top part). 



The fact that the behavior of the number variance evaluated from the traffic data coincides with the results 
obtained for interacting one-dimensional gases strengthen the idea to apply the equilibrium statistical physics 
for describing the local properties of the traffic flow. We take the advantage of this approach in a following 
thermodynamical traffic model. 

Consider N identical particles (cars) on a circle of the circumference N exposed to the thermal bath with 
inverse temperature (3. Let x\ (i = 1 . . . N) denote the position of the z-th particle and put x/v+i = %\ + N, 
for convenience. The particle interaction is described by a potential (see Ref. (3l) 

JV 

C/oc J^rf 1 , (2) 
i=l 

where rj = \xi + i — Xi\ is the distance between the neighboring particles. The nearest-neighbor interaction is 
chosen with the respect to the realistic behavior of a car-driver in the traffic stream. As published in Ref. ||5], 
the heat bath drives the model into the thermal equilibrium and the probability density Pp{r) for gap r among 
the neighboring particles corresponds to the function (Q]). 

According to HI, the number variance A n (L) of an one-dimensional gas in thermal equilibrium can be 
exactly determined from its spacing distribution Pp{r). For the clearance distribution (Q} we obtain (see iTTOlO 

A n (L)« X L + 7 (L-oo), (3) 

i.e. a linear function with a slope 

3 - 1QB/3 - 

x ~ 1 + — ; ^=^2 — 

(3 + 4^B^) 
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Figure 3 : The slope \ an d the inverse temperature (3 as functions of the traffic density 
q. The squares on the upper subplot display the slope of the number variance A n (L) 
(see Fig. 2), separately analyzed for various traffic densities. The lower subplot visual- 
izes the fitted values of the inverse temperature f3, for which the exact form of number 
variance A n (L) = x{0) L + lifi) corresponds to the number variance obtained from 
the traffic data. The continuous curves represent a polynomial approximations of the 
relevant data. 



and 

32S/3 (21 + SB 13 + 24/B/7) - 48/B^ - 63 
'■y 

245 + 

which depend on the inverse temperature j3 only. Above relations represent a large approximations whereas 
two phenomenological formulae 

1 

X ~ 2.4360 ^°- 8207 + 1 

and 



7 ~ 5.1926 /3 + 2.3929 

specify the behavior of \ an d 7 near the origin. We emphasize that, in the limiting case = 0, the value of x is 
equal to one, i.e. A n (L) = L, as expected for the independent events. The slope \ is a decreasing function of j3. 

The described properties of the function A n (L) agree with the behavior of the number variance extracted 
from the traffic data (see Fig. 2). A comparison between traffic data number variance and the formula © allows 
us to determine the empirical dependence of inverse temperature on traffic density g. The inverse temperature 
reflects the microscopic status in which the individual vehicular interactions influence the traffic. Conversely, 
in the macroscopic approach, traffic is treated as a continuum and modelled by aggregated, fluid-like quantities, 
such as density and flow (see Q). Its most prominent result is the dependence of the traffic flow on the traffic 
density - the fundamental diagram. 

It is clear that the macroscopic traffic characteristics are determined by its microscopic status. Conse- 
quently there should be a relation between the behavior of the fundamental diagram and that of the inverse 
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Figure 4: Traffic flow J(g) and inverse temperature (3(g) as functions of a traf- 
fic density g. Plus signs display a traffic flow in thousands of vehicles per hour and 
squares correspond to inverse temperature of the traffic gas. The results of a polyno- 
mial curve-fitting are visualized by the continuous curves. The dashed lines represent 
a linear approximations of the relevant data near the origin. 



temperature (3. On the Figure 4 we display the behavior of the inverse temperature (3 simultaneously with 
the fundamental diagram. The both curves show a virtually linear increase in the region of a free traffic 
(up to q ps 10 veh/km/lane). The inverse temperature j3 then displays a plateau for the densities up to 
18 veh/km/lane while the flow continues to increase. A detailed inspection uncovers, however, that the in- 
crease of the traffic flow ceases to be linear and becomes concave at that region. So the flow is reduced with 
respect to the outcome expected for a linear behavior - a manifestation of the onset of the phenomenon that 
finally leads to a congested traffic. For larger densities the temperature (3 increases up to g <; 32 veh/km/lane. 
The center of this interval is localized at g ps 25 - a critical point of the fundamental diagram at which the flow 
starts to decrease. This behavior of the inverse temperature is understandable and imposed by the fact that the 
drivers, moving quite fast in a relatively dense traffic flow, have to synchronize their driving with the preceding 
car (a strong interaction) and are therefore under a considerable psychological pressure. After the transition 
from the free to a congested traffic regime (between 40 and 50 veh/km/lane), the synchronization continues 
to decline because of the decrease in the mean velocity leading to decreasing 0. Finally - for densities related 
to the congested traffic - the inverse temperature increases while the flow remains constant. The comparison 
between the traffic flow and the inverse temperature is even more illustrative when the changes of the flow are 
taken into account. Therefore we evaluate the derivative of the flow 

J' = — 
dg 

and plot the result on the Figure 5. The behavior of the inverse temperature (3 can be understood as a quantita- 
tive description of the alertness required by the drivers in a given situation. 

The dependence of (3 on the density g can be obtained also using the measured clearance distribution and 
comparing it with the formula ([TJ). It leads to the same results as (3 obtained from the number variance A n (L). It 
is known (see (TJ) that for one-dimensional gases in thermal equilibrium the function A n (L) can be determined 
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Figure 5: Inverse temperature (3(g) and scaled derivative of the flow J' = ^ as 
functions of a traffic density g. Squares correspond to the inverse temperature of traf- 
fic sample while plus signs display the first derivative of the flow (scaled for better 
visualization). The continuous curves represent a relevant polynomial approximations. 



from the knowledge of the spacing distribution Pp(r). So the fact that obtaining j3 by virtue of the number vari- 
ance A n (L) and the spacing distribution Pp(r) leads to the same results supports the view that locally the traffic 
can be described by instruments of equilibrium statistical physics. 

In summary, we have investigated the statistical variances of the single- vehicle data from the Dutch freeway 
A9. Particularly, we have separately analyzed the number variance in eighty five equidistant density-subregions 
and found a linear dependence in each of them. Using the thermodynamical model presented originally in 
Ref.[3], we have found an excellent agreement between the number variance of particles in thermal-equilibrium 
and that of the traffic data. It was demonstrated that the inverse temperature of the traffic sample, indicating 
the degree of alertness of the drivers, shows an increase at both the low and high densities. In the intermediate 
region, where the free flow regime converts to the congested traffic, it displays more complex behavior. 

The presented results support the possibility for applying the equilibrium statistical physics to the traffic 
systems. It confirms also the hypothesis that short-ranged forwardly -directed power-law potential (0 is a good 
choice for describing the fundamental interaction among the vehicles. 
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